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SYSTEM AND METHOD FOR EFFICIENT SIMULATION OF REFLECTOMETRY 
RESPONSE FROM TWO-DIMENSIONAL GRATING STRUCTURES 

TECHNICAL FIELD 



[0001] The present invention relates to the simulation of the reflectometry response from 
grating profiles and more particularly to an efficient method for accurately simulating the 
integrated reflectometry response from two-dimensional grating structures using a few points. 



BACKGROUND ART 



[0002] Spectroscopic reflectometry and ellipsometry have been a mainstay for thin film 
U metrology for many years. Recently, spectroscopic reflectometry and ellipsometry have been 

.: ;;s. 

O applied to characterizing patterned structures in integrated circuit (IC) processing by directing a 
: beam of light on the patterned structures at a certain angle of incidence and measuring the 

:: p 

spectra of the reflected light. However, due to the difficulty of rigorous simulation responses 

i? 

y from patterned structures, empirical methods, such as neural networks (NN) or principle 

t*\ 
■■■ iip 

f || component analysis (PC A), are used to build up the relation between the patterned structure 

=: U 

Q parameters (e.g., width of a structure (CD), grating height, sidewall angle, and other profile 

: : 

si* 

parameters) and reflected spectra. 



[0003] Figure 1 illustrates a typical reflectometry configuration for characterizing patterned 
structures located on a wafer. A broadband light beam 105 travels through an optical system 1 10 
characterized by a numerical aperture. A lens 115 focuses the broadband light beam 105 into a 
spot characterized by a spot size onto a patterned structure located on wafer 120 (note that 
throughout this specification a lens is referred to more genetically as an aperture). The light 
reflected off of the patterned structure located on wafer 120 is then collected by lens 115 and 
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transmitted to a spectrometer 125 through optical system 110. Using other metrology tools, such 
as a scanning electron microscope (SEM), an atomic force microscope (AFM), etc., the 
parameters of the patterned structure can be measured off-line. Next, neural networks or 
principal component analysis can be used to build up (train) a non-linear relation between the 
reflected spectra and the parameters of the patterned structure. However, the relation between 
the reflected spectra and the parameters of the patterned structure is only valid when the values 
of the parameters are within the range of parameters used for training. Furthermore, long 
turnaround times imposed by off-line metrology also result in lower yields, slower learning 
curves, and higher costs for new processes and products. 

[0004] A multi-point rigorous simulation method can be used to avoid doing experiments to 
build up empirical relations between the reflected spectra and the parameters of patterned 
structures. Ideally, the numerical aperture should be very small (e.g. less than 0.01) so that the 
reflectometry response can be simulated using normal incidence. However, if the lens 1 15 is too 
small, then the light throughput is low, thus weakening the intensity of the light beam 105 
incident on a two-dimensional grating structure. The weaker the intensity of the light beam 105, 
the longer it takes to collect enough simulation data to achieve stability. 

[0005] Starting from Maxwell's equations, the response of light reflected from patterned 
structures can be simulated rigorously using numerical methods. However, in an actual 
reflectometry system, to have an acceptable spot size (less than 100 |tim) and throughput (less 
than 1 second), the numerical aperture (NA) is fairly large (about 0.05 or larger). Thus, the 
reflectance is actually the integrated response from multiple beams of light reflected off the 
patterned structures, where the incidence angle of each light beam is close to zero degrees. 
Figures 2a and 2b show an example of the multi-point rigorous simulation method. Figure 2a 



shows light passing through numerous (e.g., 20 to 30) points 205 across an aperture 210 creating 
numerous light beams 215. The aperture 210 focuses each light beam 215 onto the two- 
dimensional grating structure 220 at an angle close to zero degrees. Simulating the numerous 
light beams 215 focused onto the two-dimensional grating structure 220 yields the response 
distribution 225 for light beams 215 as shown in Figure 2b. The integrated response for this 
wavelength very closely approximates the actual reflectance and can be obtained by adding up 
each of the responses shown in Figure 2b. However, the simulation takes a long time. 

[0006] A single point rigorous simulation method can also be used to avoid doing 
experiments to build up empirical relations between the reflected spectra and the parameters of 
the patterned structures. Figures 3a and 3b show an example of the single point rigorous 
simulation method. Figure 3a shows a simulation of light passing through a single point O 305 
located at the center of an aperture 310 creating a single light beam 315. The light beam 315 is 
incident on the two-dimensional grating structure 320 at an angle of zero degrees. Fig. 3b shows 
the simulated reflectance response 325 of the light passing through the single point O 305. 
Although, the reflectance response 325 can be obtained very quickly, it is not very accurate since 
it only represents the reflectance response 325 from light passing through the single point O 305. 
The reflectance response 325 does not take into account the reflectance response from light 
passing through other points across the aperture 310, and thus cannot represent the overall 
reflectance response of the aperture 310. 

[0007] Figure 4 is a simulation which compares the accuracy of the response from light 
reflected off of a two-dimensional grating structure using the multi-point rigorous simulation 
method and the single point simulation method. As shown in Figure 4, the simulated response 
410 obtained using the single point rigorous simulation method differs greatly from the response 



405 obtained using the multi-point rigorous simulation method. The multi-point rigorous 
simulation method closely approximates the actual reflectance response. Thus, there is a desire 
for a method of simulating the reflectance response of two-dimensional grating structures that is 
both accurate and not time consuming. 

SUMMARY OF INVENTION 

[0008] The method in accordance with embodiments of the present invention relates to a 
method for efficient simulation of reflectometry response from two-dimensional grating 
structures. 

[0009] In one embodiment, the light intensity distribution across the aperture is uniform. A 
first and a second point are determined within an aperture located in an optical system. Next, the 
reflectance response of light incident at the first point and the second point are simulated. The 
approximated integrated reflectance response of the aperture is then determined based on the 
reflectance response at the first point and the second point and determined characteristics of the 
optical system. 

[0010] In another embodiment, the light intensity distribution across the aperture is not 
uniform. A first and a second point are determined within an aperture located in an optical 
system. Next, the reflectance response of light incident at the first point and the second point are 
simulated. The approximated integrated reflectance response of the aperture is then determined 
based on the reflectance response at the first point and the second point and determined 
characteristics of the optical system. 



BRIEF DESCRIPTION OF DRAWINGS 

[0011] FIG. 1 shows a typical reflectometry configuration for characterizing patterned 
structures located on a wafer. 

[0012] FIG. 2a illustrates the prior art multi-point rigorous simulation method for simulating 
the response from light reflected off of a two-dimensional grating structure. 

[0013] FIG. 2b shows a plot of the reflectance response obtained using the multi-point 
rigorous simulation method illustrated in Fig 2a. 

[0014] FIG. 3a illustrates the prior art single point rigorous simulation method for simulating 
the response from a single light beam reflected off of a two-dimensional grating structure. FIG. 
3b shows a plot of the reflectance response obtained using the single point rigorous simulation 
method illustrated in Fig 3a. 

[0015] FIG. 4 shows a simulation comparing the accuracy of the simulated response from 
light reflected off of a two-dimensional grating structure using the multi-point rigorous 
simulation method and the single point rigorous method. 

[0016] FIG. 5a shows a lens focusing a light beam passing through point x z onto a two- 
dimensional grating structure, where the angle of the light beam incident on the two-dimensional 
grating structure is 6 t . FIG. 5b shows an example of a 2-angle quadrature method for simulating 
the integrated response from light reflected off of a two-dimensional grating structure. 

[0017] FIG. 6 shows a three dimensional view of an aperture having two incidence angles 
along the line GG' which is parallel to a two-dimensional grating structure. 



[0018] FIG. 7 shows a simulation comparing the reflectance response of light passing 
through two points, sj and S2, along a line parallel to a two-dimensional grating structure. 

[0019] FIG. 8 shows a simulation comparing the simulated response from light reflected off 
of a two-dimensional grating structure using the multi-point rigorous simulation method and an 
approximation of the simulated response curve using a 6th order even polynomial function. 

[0020] FIG. 9 shows a flow chart overview of a process to derive the Wangle quadrature for 
numerical aperture integration for a two-dimensional grating structure where the light intensity 
distribution across the aperture is uniform. 

[0021] FIG. 10 shows a simulation comparing the accuracy of the simulated response from 
light reflected off of a two-dimensional grating structure using the multi-point rigorous 
simulation method, the single point rigorous method, and the 2-angle quadrature method. 

[0022] FIG. 1 1 shows an example of the 2-angle quadrature method for simulating the 
integrated response from light reflected off of a two-dimensional grating structure, where the 
light intensity distribution across the aperture is not uniform. 

[0023] FIG. 12 shows a flow chart overview of a process to derive the £- angle quadrature for 
numerical aperture integration for a two-dimensional grating structure where the light intensity 
distribution across the aperture is not uniform. 

DETAILED DESCRIPTION 

[0024] The present invention provides a method for efficient simulation of reflectometry 
response from two-dimensional grating structures. In the following description, numerous 
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details are set forth in order to enable a thorough understanding of the present invention. 
However, it will be understood by those of ordinary skill in the art that these specific details are 
not required in order to practice the invention. Further, well-known elements, devices, process 
steps and the like are not set forth in detail in order to avoid obscuring the present invention. 

[0025] Although the invention has been described in conjunction with particular 
embodiments, it will be appreciated that various modifications and alterations may be made by 
those skilled in the art without departing from the spirit and scope of the invention. The 
invention is not to be limited by the foregoing illustrative details, but rather is to be defined by 
the appended claims. 

[0026] As stated above, there are a number of ways to simulate the reflectometry response of 
the two dimensional grating structure. However, each of the prior art methods for simulating the 
reflectometry response of the two dimensional grating structure are either too costly, too time 
consuming, not accurate enough, or a combination thereof. Thus, there is a desire for an efficient 
method for simulating the reflectometry response of two-dimensional grating structures that 
overcomes the problems associated with the prior art. The present invention discloses a method 
for deriving the Wangle quadrature for numerical aperture integration for a two-dimensional 
grating structure, where k is the number of simulation angles used to approximate the 
reflectometry response of the two-dimensional grating structure. 

[0027] Figure 5a shows an aperture 505 focusing a light beam 508 passing through point x t 
onto a two-dimensional grating structure 510. The angle of the light beam 508 incident on the 
two-dimensional grating structure 510 is 0 f . The numerical aperture of the optical system is 
defined as «(sin 0 m ), where n is the refractive index of the medium and 6 m is the half-angle of the 



cone 515 formed as light is focused through the aperture 505 onto the two-dimensional grating 
structure 510. The test grating area for many integrated circuits is generally only 100 fxm or less 
in size. Therefore, the spot size of the light beam 508 incident on the two-dimensional grating 
structure 510 located within the test grating area must also be very small. If the spot size of the 
light beam 508 is larger than the test grating area, light may reflect off other features located on 
the integrated circuit, thus distorting the desired reflectance response of the two-dimensional 
grating structure 510. Generally, the value of the half-angle d m is determined by the 
measurement tools used to measure the two-dimensional grating structures 510 and can be 
obtained from the measurement tool manufacturer. 

[0028] To simplify the math used to derive the fc-angle quadrature for numerical aperture 
integration for a two-dimensional grating structure 510 in a first embodiment of the present 
invention, some assumptions and premises are made. First, an assumption is made that the 
intensity of the light incident on the aperture 505 is uniformly distributed. The information as to 
whether or not the light distribution is actually uniformly distributed across the aperture 505 can 
be obtained from the manufacturer of the metrology tool. Note, however, that in a second 
embodiment of the present invention, if the actual light intensity distribution across the aperture 
505 is not uniform, but can be measured, the measured light intensity distribution can be used as 
weight in determining the Wangle quadrature for numerical aperture integration for a two- 
dimensional grating structure. The second embodiment of the present invention is explained 
further below. 

[0029] To further simplify the math used to derive the Wangle quadrature, a first premise is 
made that if the two-dimensional grating structure 510 is substantially parallel (which includes 
parallel) to the line GG 5 525 located on the aperture 505 as shown in FIG. 5b, then the 
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reflectance response of each light beam 508 passing through any point along the line GG' 525 
and onto the two-dimensional grating structure 510 is approximately the same. This premise can 
be verified using rigorous simulation on a computer (see discussion of Figures 6 - 7 below). 



[0030] Fig. 6 shows a three-dimensional view of the aperture 505 having a line GG' 525 
parallel to the two-dimensional grating structure 510 (shown in Fig. 5a) and a horizontal line LL' 
530 across the center of the aperture 505 and perpendicular to the line GG' 525. Two points si 
and S2 are also shown along the line GG' 525, where S2 is located at the intersection of the lines 
LI/ 530 and GG' 525. The incident angle of the unpolarized light beam 605 passing through the 
point sj is four degrees and the incident angle of the unpolarized light beam 610 passing through 
the point S2 is two degrees. To verify the first premise, the reflectance response of the light beam 

; ; us. 

1* 605 passing through the point sj and the light beam 610 passing through the point s 2 are 
i,ft simulated using the multi-point rigorous simulation method and shown in Fig. 7. 

i 

M [0031] Figure 7 shows that the reflectance response 705 of the light passing through point sj 

□ 

HI! and the reflectance response 710 of light passing through point S2 located on along the line GG' 

US 

•lit* 

?" 1 525 are almost the same, thus validating the first premise stated above. The reflectance 
responses for light beams passing through other points along the line GG' 525 can also be 
simulated and compared to further validate this premise. Therefore, the reflectance response of 
light passing through only one point on the line GG' 525 approximates the reflectance response 
for the entire line GG' 525 if the line GG' 525 is parallel to the two-dimensional grating 
structure 510. This is true for any line within the aperture 505 that is parallel to the two- 
dimensional grating structure 510 and if needed, can be verified using rigorous simulation. 
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[0032] In one embodiment, the point chosen to simulate the reflectance response for each 
line within the aperture 505 parallel to the two-dimensional grating structure 510 is located at the 
center of each respective line. For example, the preferred point for approximating the reflectance 
response for the line GG' 525 is the point s 2 . Note, however, that each point located at the center 
of each line within the aperture 505 parallel to the two-dimensional grating structure 510 will 
have a different reflectance measurement. Thus, all that is required to approximate the 
reflectance of the entire aperture 505 is to simulate the reflectance response of light passing 
through every point along the line LL' 530. 

[0033] To even further simplify the math used to derive the Wangle quadrature, a second 
premise is made that the reflectance response R(x) for any wavelength from the two-dimensional 
grating structure 510 along the direction of the line LL' 530 can be approximated by using a set 
of orthonormal functions, such as polynomials, Fourier series, or any other orthogonal and 
normal function sets. For example, since the reflectance response 525 is continuous and 
symmetrical about the point O 503, a (4£-2) th order even polynomial function can be used. If 
k=2, then R(x) is a 6 th order polynomial even function as shown in equation (1) below: 



R(x) = a 0 +a. 



"V" -v- \ / \ 



X 



\ X mJ 



x 



+ a. 



[0034] Where a t (i = 0,1,2. ..) depends on wavelength and the two-dimensional grating 
structure 510. The number of orders needed depends on how fast the reflectance R(x) changes 
along the line LL' 530. For a numerical aperture of a typical reflectometer (e.g., approximately 
0.1 or less) and most two-dimensional grating structures, this condition is satisfied for k = 2 or 3. 
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[0035] This second premise can be verified by simulating the response from multiple light 
beams (e.g., 10) reflected off of the two-dimensional grating structure 510 using the prior art 
multi-point rigorous simulation method on a computer, where the light beams focused onto the 
two-dimensional grating structure 510 pass through multiple points (e.g., 10) located along the 
line LL' 530. Once the reflectance response from the multiple points are simulated, the next step 
is to approximate the resulting curve of the reflectance response with a (4£-2) th order even 
polynomial function as illustrated in Figure 8 for£= 2. As stated above, for a numerical 
aperture of a typical reflectometer (approximately 0.1 or less) and most two-dimensional grating 
structures, this second premise is satisfied for k = 2 or 3. However, it should be noted that k is 
preferably equal to a value sufficient to approximate the polynomial for the resulting wavelength. 
Thus, if the second premise is not satisfied for k = 2 or 3, higher values ofk can be used. 

[0036] Fig. 9 shows a flow chart overview 900 of a process to derive the Wangle quadrature 
for numerical aperture integration for a two-dimensional grating structure. The first step 905 in 
deriving of the Wangle quadrature for numerical aperture integration for a two-dimensional 
grating structure is to obtain the numerical aperture of the optical system and the light intensity 
distribution across the aperture 505 from the manufacturer of the metrology tool. Note that the 
numerical aperture and the light intensity distribution across the aperture 505 can also be 
measured. For purposes of this embodiment, it is assumed that the light distribution within the 
aperture 505 is uniform, the numerical aperture of the optical system is approximately 0.1 or less, 
and k = 2, i.e., the 2-angle quadrature is used to approximate the integrated reflectance response. 
However, in other embodiments, the numerical aperture can be greater than 0.1 and for other 
values of k, the derivation can be generalized accordingly. 
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[0037] The next step 910 is to determine two optimal points, x } and x 2 , along the line LL' 
530 from which to approximate the integrated reflectance response of the aperture 505. The 
weight of the reflectance response from the two light beams passing through optimal points xj 
and x 2 is determined to approximate the reflectance response of the entire aperture 505. The 
weights for the light intensity distribution at optimal points xj and x 2 are wj and w 2 , respectively. 
An equation using the weighted average of reflectance responses from points xi andx 2 to 
approximate the integrated reflectance across the aperture 505 is shown in equation (2) below. 

\R(x)*^xl-x 2 •4*dx = [w l R(x 1 ) + w 2 R(x 2 )]*itx 2 m , (2) 

0 

[0038] where R(x) is the reflectance at any incidence angle along the line LL' 530, 

4 x l -x 2 & is the weight of the area at position x along the line LL' 530, and x m is the length of 

the line OL'. The value of x m can be calculated from the value of numerical aperture. 

[0039] Note that making the assumption that the intensity of the light is uniformly distributed 
across the aperture 505 simplifies equation (2) by eliminating an additional term specifying the 
weight due to any non-uniform distribution of light across the aperture 505. Eliminating the 
additional term specifying the weight due to any non-uniform distribution of light across the 
aperture 505 ensures a closed-form analytical solution of the integral on the left hand side of 
equation (2). However, it should be noted that there can also be a closed-form analytical solution 
of the integral on the left hand side of equation (2) in some cases where the light distribution is 
not uniform. 

[004.0] In step 9 1 5, the first premise is made that response of light passing through only one 
point on any line within the aperture 505 that is parallel to the two-dimensional grating structure 
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510 is sufficient to approximate the reflectance response for that line within the aperture 505 if 
that line is parallel to the two-dimensional grating structure 510. Again, this first premise can be 
verified using rigorous simulation on a computer as described above. In step 920, the second 
premise is made that that the reflectance response R(x) for any wavelength from the two- 
dimensional grating structure 510 along the direction of the line LL' 530 can be approximated by 
a set of orthonormal functions, such as polynomials, Fourier series, or any other orthogonal and 
normal function sets. For example, since the reflectance response 525 is continuous and 
symmetrical about the point O 503, a {Ak-lf order even polynomial function can be used. If 
k=2, then R(x) is a 6 th order polynomial even function as shown in equation (1) above. Again, 
the second premise can be verified by simulating the response from multiple light beams (e.g., 
10) reflected off of the two-dimensional grating structure 510 using the prior art multi-point 
rigorous simulation method on a computer as described above. 

[0041 ] In step 925, substitute R(x) from equation (1) into R(x) of equation (2) and then 
perform the integration which yields equation (3) below: 



( w i + ^2 - lk + ^Wi + -£)« 2 + 



]«4 +[ ^i 6 + ^ 2 6 -^-|fl 6 =0 (3) 



64, 



where *,-=—, (i = 1,2) . Since equation (3) should be valid for any value of a, (i = 0, 1 ,2, . . .), 
their coefficients should be equal to 0 as shown below. 



w l +w 2 -1 = 0 (4) 



wtf+w 2 t 2 2 ~ = 0 (5) 
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vV. 4 + w 2 f 2 4 -| = 0 (6) 



*\ti+vrf~ = 0 (7) 

04 



[0042] In step 930, solving equations (4) through (7) yields one solution set for w u w 2 , and f, 
and t 2 having real physical meaning: 

w 1 =0.723607 
w 2 =0.276393 
t x =0.309017 
t 2 =0.809017 

[0043] However, note that solving equations (4) through (7) may yield multiple solution sets 
for w„ w 2 , and t x and t 2 having real physical meaning. Considering the relation between t t and 
incidence angle 0, yields: 



' ~ tan0„ ' (8) 



therefore, 



8,. = arctan{r. • tan[arcsin(iV^)]} (9). 
[0044] For example, if NA = 0.07, then 0, = 1 .2422, B 2 = 3 .2492. However, NA could be any 
value sufficient to make the previously stated assumptions and premises valid. 

[0045] Once the values of d,, 6 2 , wj and w 2 are calculated for a particular numerical aperture, 
the values wj, w 2 , 9 h and 6 2 can then be used to rigorously simulate the reflectance response of 
the light incident at points x u and x 2 using a computer in Step 935. The reflectance response 
values R(x 7 ) and R(x 2 ) are then obtained from the resulting rigorous simulation. Finally the 
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values for R(x } ), R(x 2 ), wj and w 2 are used to calculate the weighted average to represent the 
integrated response from the aperture 505 in step 940. 

WrRixJ + WiRixJ (10) 
[0046] Fig. 10 shows a simulation comparing the results of the simulated reflectance 
response across the aperture 505 from a two-dimensional grating structure using the prior art 
multi-point rigorous method (accurate simulation) 1005, the prior art single-point rigorous 
method (0 degree incidence) 1010, and the Wangle quadrature method 1015 in accordance with 
an embodiment of the present invention where k = 2. The plot shown in Fig. 10 shows that the 
simulated reflectance response of the aperture 505 using the 2-angle quadrature method very 
closely approximates the simulated reflectance response of the aperture 505 using the prior art 
multi-point rigorous simulation method. The simulation shown in Fig. 10 also shows that the 
simulated reflectance response of the aperture 505 using the 2-angle quadrature method is 
significantly more accurate than simulated reflectance response of the aperture 505 using the 
prior art single point rigorous simulation method. Thus, the Wangle quadrature method of the 
present invention provides a significant advantage over the prior art multi-point rigorous 
simulation method by very closely approximating the simulated reflectance response of the 
aperture 505 twenty to thirty times more quickly. Similarly, the Wangle quadrature method of 
the present invention provides a significant advantage over the prior art single point rigorous 
simulation method by quickly and accurately approximating the reflectance response of the 
aperture 505. Furthermore, the Wangle quadrature method provides a significant advantage over 
the prior art in that once the values of x u x 2 , 0 l9 6 2 , w } and w 2 are calculated for a particular 
numerical aperture, the overall integrated reflectance response across the aperture 505 can be 
determined for any two-dimensional grating structure using only ^-angles. 
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[0047] Figure 1 1 illustrates a second embodiment of the present invention where the light 
distribution within the aperture 1 105 is not uniform. Although the light distribution within the 
aperture 1 105 is not uniform, the same method for deriving the Wangle quadrature for numerical 
aperture integration for a two-dimensional grating structure disclosed in the first embodiment can 
be used if the actual light intensity distribution across the aperture 1 105 can be measured. The 
measured intensity distribution of the light incident on the aperture 1 105 is used as weight in 

equation (2) in place of the weight term ^Jx 2 m -x 2 dx . In addition, numerical integration may be 
needed. 

[0048] Figure 12 illustrates a flow chart overview 1200 of a process to derive the Wangle 
quadrature for numerical aperture integration for a two-dimensional grating structure. The first 
Step 1205 is to obtain the numerical aperture of the optical system and the light intensity 
distribution across the aperture 1 105 from the manufacturer of the metrology tool Note that the 
numerical aperture and the light intensity distribution can also be measured. For purposes of this 
example, it is assumed that the light distribution within the aperture 1 105 is not uniform, the 
numerical aperture of the optical system is approximately 0.1 or less, and k = 2, i.e., the 2-angle 
quadrature is used to approximate the integrated reflectance response. However, in other 
embodiments, the numerical aperture can be greater than 0.1 and for other values of k, the 
derivation can be generalized accordingly. 

[0049] In Step 1210, the first premise, as stated above, is used to simplify the math. If a two- 
dimensional grating structure is parallel to the line GG 5 1110, then the reflectance response of 
each light beam passing through any point along the line GG' 1 1 10 and onto the two- 
dimensional grating structure is almost the same. Likewise, if the two-dimensional grating 
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structure is parallel to the line HH' 1115, then the reflectance response of each light beam 
passing through any point along the line HH' 1115 and onto the two-dimensional grating 
structure is almost the same. This first premise can be verified using rigorous simulation on a 
computer as discussed in Figures 6 - 7 above. Therefore, the reflectance response of light 
passing through only one point is sufficient to approximate the reflectance response for the entire 
line GG' 1 1 1 0 or line HH' 1115. This first premise is true for any line within the aperture 1 1 05 
that is parallel to a two-dimensional grating structure. In a preferred second embodiment, the 
point chosen to simulate the reflectance response for each line within the aperture 1 105 parallel 
to the two-dimensional grating structure is located at the center of each respective line. For 
example, the preferred point for approximating the reflective response for the line GG' 1 1 10 is 
the point xj and the preferred point for approximating the reflective response of the line HH' 
1 1 1 5 is the point x 2 . Note, however, that each point located at the center of each line within the 
aperture 1 105 parallel to the two-dimensional grating structure will have a different reflectance 
measurement. 

[0050] In Step 1215, the second premise, as stated above, is used again to even further 
simplify the math used to derive the Wangle quadrature. The reflectance response R(x) for any 
wavelength from the two-dimensional grating structure along the direction of the line LL' 1 125 
can be approximated by a set of orthonormal functions, such as polynomials, Fourier series, or 
any other orthogonal and normal function sets. In this particular embodiment, a (4£-2) th order 
even polynomial function can be used since the reflectance response is continuous and 
symmetrical about the point O 503. For example, if £=2, then R(;c) is a 6 th order polynomial even 
function as shown in equation (1 1) below: 
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R(x) = a 0 +a 7 
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(11) 



[0051] Where a, (i = 0,1,2. . .) depends on wavelength and the two-dimensional grating 
structure. The number of orders needed depends on how fast the reflectance R(x) changes along 
the line LL' 1 125. For a numerical aperture of a typical reflectometer (e.g., approximately 0.1 or 
less) and most two-dimensional grating structures, this condition is satisfied for £ = 2 or 3. This 
second premise can be verified as discussed above and shown in Figure 8. 

[0052J The next Step 1220 is to determine two optimal points, xi and x 2 , along the line LL' 
1 125 from which to approximate the integrated reflectance response of the aperture 1 105. The 
weight of the reflectance from the two light beams passing through the optimal incidence angles 
xj and x 2 is also be determined to approximate the reflectance response of the entire aperture 
1105. 

[0053] As stated above, the light distribution incident on the aperture 1 105, is not uniform, 
but can be measured. The non-uniform distribution of light on aperture 1 105 can be specified 
with a function l(x,y) or l(r,Q ), where r and 0 are the radial and the angular coordinates, 

respectively. Thus, the weight term <Jx 2 m -x 2 dx on the left hand side of the integral in equation 
(2) can be replaced with a new weight term w(x)dx representing the measured intensity 
distribution of the light incident on the aperture 1 105. 

[0054] The weight term w(x)dx representing the measured intensity distribution of the light 
incident on the aperture 1 105 is equal to the amount of the light passing through line GG' 1110 
having a width dx plus the amount of light passing through the line HH' 1 1 1 5 having a width dx. 
Note that line GG' 1110 located at position x } and line HH' 1115 located at position x 2 are both 
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symmetric about the center 1 120 of the aperture 1 105. Therefore, weight term w(x)dx can be 
written as 

w(x)dx=Tix 2 m dx \l{x,y)dy, (12) 

where the light distribution function l{x,y) satisfies the normalization condition 

\\l{x,y)dxdy = \. (13) 

aperture 

[0055] For example, if the light distribution function is a circularly symmetric function: 



l(x,y) = l(r,Q) = ± 



2x1 



, (14) 



where the normalization factor kisK = -nx 2 m , then due to the circular symmetry of the light 



distribution function, the weight term is: 

4*1-* 



w(x)dx = 4nx 2 m dx jl(x, y)dy = dx^x 2 m -x 



6 3x 2 



(15) 



m J 



[0056] Replacing the weight term ^x 2 m -x 2 dx on the left hand side of equation (2) with the 
weight term in equation (15) above yields: 



jR(xH*)dx=jnx 2 m jR( X yx 2 m -x< 



2 \ 



5_l_x_ 
6 3^ y 



dx. (16) 

where R(x) is the reflectance for the light beam passing through point Jtalong the line LL' 1 125, 



w{x)dx = ^dx^x 2 m -x' 
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[0057] is the weight of the area at position x along the line LL' 1 1 25, and x m is the length of 
the line OL' 1125. 

[0058] In Step 1225, substitute R(x) from equation (11) into R(x) of equation (16) and then 
performing the integration yields equation (18) below: 



(w, + w 2 -l)a 0 +j\v,;, 2 + w 2 t 2 2 ~a 2 + ^w,t} 



w,f, s +w 2 r 2 6 ~k=0 (18) 



where /, = — l -, (i = 1,2) . Since the above equation should be valid for any value of a, 
(i = 0,1,2,.. .), their coefficients should be equal to 0 as shown below. 

w,+w 2 -l = 0 (19) 



wtf+w 2 t 2 2 ~ = 0 (20) 



^:+w 4 2-—=o (2i) 



^1+^2-^ = 0- (22) 



[0059] In Step 1230, solving equations (19) through (22) yields one solution set for w„ w 2 , 
and t\ and t 2 having real physical meaning: 



w, =0.247517 
w 2 = 0.752483 
t x =0.296618 
t 2 = 0.793932 

[0060] However, note that solving equations (19) through (22) may yield multiple solution 
sets for w„ w 2 , and r, and t 2 having real physical meaning. Considering the relation between t, 
and incidence angle 0, yields: 
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tanO. 

(23) 

tan9 m 



therefore, 



0 ; = arctanfo • tan[arcsin(A^)]} . (24) 
[0061] For example, if NA = 0.07, then 6 } = 1.1924, 6 2 = 3.1888. However, AM could be any 
value sufficient to make the previously stated assumptions and premises valid. 

[0062] Once the values of <9 ; , 6 29 wj and w 2 are calculated for a particular numerical aperture, 
the values wj 9 w 2 , 6 J9 and 6 2 can then be used rigorously simulate the reflectance responses of 
light beams passing through points x h and x 2 using a computer in Step 1235. The reflectance 
response values R(x } ) and R(x 2 ) are then obtained from the resulting simulation. Finally the 
values for R(x } \ R(x 2 ), w } and w 2 are used to calculate the weighted average to represent the 
integrated response from the aperture 505 in step 1240. 

w^J + vv^) (25) 
[0063] Thus, the Wangle quadrature method of the present invention provides a significant 
advantage over the prior art multi-point rigorous simulation method by very closely 
approximating the simulated reflectance response of the aperture 1 105 twenty to thirty times 
more quickly. Similarly, the Wangle quadrature method of the present invention provides a 
significant advantage over the prior art single point rigorous simulation method by quickly and 
accurately approximating the reflectance response of the aperture 1 105. Furthermore, the k- 
angle quadrature method provides a significant advantage the over the prior art in that once the 
values of xj 9 x 2t 6 h 6 2y wj and w 2 are calculated for a particular numerical aperture, the overall 
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integrated reflectance response across the aperture 1 105 can be determined for any two- 
dimensional grating structure using only ^-angles. 



fill 
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